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ABSTRACT 

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880 fim 
and the Keck Adaptive Optics (AO) system at -Ks-band of a gravitationally lensed sub-millimeter 
galaxy (SMG) at z = 4.243 discovered in the Herschel Astrophysical Terahertz Large Area Survey. 
The SMA data (angular resolution » 0'.'6) resolve the dust emission into multiple lensed images, 
while the Keck AO its-band data (angular resolution « O'.'l) resolve the lens into a pair of galaxies 
separated by 0'.'3. We present an optical spectrum of the foreground lens obtained with the Gemini- 
South telescope that provides a lens redshift of zi cns = 0.595 ± 0.005. We develop and apply a new 
lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of 
fi = 4.1±0.2 and has an intrinsic infrared (IR) luminosity of -Lir = (2.1±0.2) x 10 13 L & . We measure a 
half-light radius of the background source of r s = 4.4 ±0.5 kpc which implies an IR luminosity surface 
density of Sir = (3.4 ± 0.9) x 10 11 i© kpc 2 , a value that is typical of z > 2 SMGs but significantly 
lower than IR luminous galaxies at z ~ 0. The two lens galaxies are compact (n ons w 0.9 kpc) 
early-types with Einstein radii of #ei = 0.57 ± 0.01 and 6*E2 = 0.40 ± 0.01 that imply masses of 
Miensi = (7.4 ± 0.5) x 10 10 M Q and Mi ens2 = (3.7 ± 0.3) x 10 10 M . The two lensing galaxies are 
likely about to undergo a dissipationless merger, and the mass and size of the resultant system should 
be similar to other early-type galaxies at z ~ 0.6. This work highlights the importance of high spatial 
resolution imaging in developing models of strongly lensed galaxies discovered by Herschel. 
Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: high-redshift 
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1. INTRODUCTION 

It has been known for over a decade that the star- 
formation rate density i n the Universe peaked around 
redshifts z = 1 - 3 (e.g., iMadau et aLlll996t iLillv et al.l 
1996). More recently, the advent of bolometer ar- 
rays in the sub-millimeter (sub-mm) as well as the 
Spitzer Space Telescope, have established that the con- 
tribution of ultra-luminous infrared galaxies (ULIRGs) 
to the star-formation rate density in the Universe 
rises sharply w i th redshift out to z ~ 2 (e.g., 
Blain et all 119991; IChapman et al.l 120051; iLe Floc'h et al.l 
20051: iMurphv et all 1201 U iMagnelli et al.l 1201 ![) " Ah 



thoug h ULIRGs in the local Universe are known to be 
rare (|Soifer et al.l 119861 ) and have long been thought 
to arise from a major merger of two gas-rich disk 
galaxies (e.g.. JArmus et al.l 119871: IMurphv et al.l 119961 : 
iClements et aLf ll996; Bus house et al.ll2002D . their nature 
and role in galaxy evolution at high redshift is not yet 
well- understood. 

The primary obstacle to studying ULIRGs at high- 
redshift has been one of identification (caused in 
large part by faintness at optical wavelengths). Sur- 
veys to identify ULIRGs have either been limited to 
small areas on the sky (e.g., the Sub-mm Common 
User Bolometer A rray Half Degree Survey, SHADES; 
iCoppin et al.1120 06) . low-spatial resolution imaging (e.g., 
the Ba lloon-borne Large A perture Submillimeter Tele- 
scope; iPascale et all l2008h or are sensitive to mid- 
infrared (mid-IR) radiation which is far from the peak 
of the spectral energy distr ibution (SED) of the ULIRG 
(e.g., lLonsdale et al l [20031 ). Each of these techniques 
produce samples of objects that are sufficiently faint 
at far-IR and sub-mm wavelengths that follow-up ob- 
servations have been time-consuming and therefore lim- 
ited to a modest number o f objects, both in term s of 
determining redshifts (e.g.. IChapman et al.l I2005T) and 
measuring important quantities such as accurate posi- 
tions (e.g. JDannerbauer et al.ll2002tlYoung cr 2007) , mor- 
pholo g ies fe.g.. IBussmann et al.ll2009at iSwinban k et al.l 
120101 : IBussmann et al.l 1201 ID. and gas and dust 
masses fe.g iGreye et al.l 12005] ITacconi et al.l [20061: 

ilr 



Coppin et al.ll2008l: ITacconi et al 



2009b : iMichalowski et al.l 12010 ; 



120081: IBussmann et al l 
Kovacs et al.l 120101: 



Ivison et al.ll201ll : iRiechers et alJ[201ll ). 

This situation is now being remedied following the 
launch of the Herschel Space Observatory {Herschel). 
With a large array of sensitive detectors at 70 /im, 
100 /im, 160 /im, 250 /im, 350 /im, and 500 /im, the Pho- 
todetector Array Ca mera and Spectrometer (PACS; 
Poglitsch et al.l 120101 1 and Spectral and Ph otometric 
Imaging Receiver (SPIRE; iGriffin et all 120101 ) on Her- 
schel are well-suited to surveying large areas of the sky 
at wavelengths that are ideal for the detection of ULIRGs 
in the redshift range z ~ 2 — 4. The widest such survey 
is known as the Herschel Astrophysical Ter ahertz Large 
Area Survey (H- ATLAS; lEales et al.l 120101 ) and covers 
550 deg 2 of sky as the largest open-time key project, 
reaching 5<r sensitivities of 130 mJy at 100 /im, 120 mJy 
at 160 fj,m, 32 mJy at 250 mm, 36 mJy at 350 /im and 
45 mJy at 500 /im (jlbar et al.ll2010L IPascale et al.ll2011fc 
IRjgbv et al.ll2T)lTI ). 

The wide area coverage of H- ATLAS makes it ideal for 
building statistically significant samples of rare galax- 



ies. One such example that has been particularly fruit- 
ful thus far is the selection of gravitationally lensed ob- 
jects. These are systems where the light from a distant 
source (in this ULIRG at z ~ 4) is deflected 

by a foreground lens (typically an early type galaxy or 
group of galaxies) in such a way that the background 
ULIRG appears to have its angular size and bright- 
ness increased. Several authors have predicted that the 
sub-mm is an efficient waveband to identify lensing sys- 
tems due to the steep number counts of galaxies se- 
lected at sub-mm and mm wave lengths (SMGs) (e.g ., 
[Blainl [1991 iPerrotta et all 12001 iNegrello et all 120071 ). 
Additionally, the fact that most SMGs lie at z > 2 
(|Chapman et al.l 12005) increases the probability that an 
interlo ping galaxy wil l lie a long the line-of-sight. Re- 
cently, Negr ello et al.l (|2010f ) have shown that a selec- 
tion at 500 /im of -FsoOjum > 100 m Jy sources within 
the 14.4 deg 2 Science Demonstration Phase field of H- 
ATLAS identifies strongly lensed s ystems, low-z spiral 
galaxies ()Serieant fc Harrisonll2005D . and higher- z active 
galactic nuclei (AGN) that are radio-bright and show 
a sync hrotron emission sp ectrum even into the SPIRE 
bands (|de Zotti et al.ll2"o"05T ) . Shallow ground-based opti- 
cal and radio imaging can be used to remove the latter 
two classes of objects, leaving only the strongly lensed 
systems. 

The H-ATLAS source catalog already extends to w 
130 deg 2 (Phase 1 catalog; Dunne et al., in prep.). This 
paper focuses on one source of particular interest drawn 
from the Phase 1 catalog: H-ATL AS J142413.9+0 22304 
(this object is denoted "ID 141" in lCox et al.ll2011l here- 
after, we refer to it as G15v2.779). SPIRE photometry 
of this source shows that it is one of the brightest de- 
tected so far in Herschel wide-field surveys and that its 
SED peaks at wavelengths greater than 500 /im, suggest- 
ing that it lies at z > 3. This source has been the target 
of significant follow-up efforts: the Plateau de Bure In- 
terferometer (PdBI) has detected millimeter (mm) and 
sub-mm CO emission lines which imply that the red- 
shift of this source is z = 4.243 ± 0.001, while data 
from the Atacama Pathfinder Experiment (APEX) have 
shown that the dominant coolin g line in this galaxy is 
[CII] emission ()Cox et al.l 1201 1|) . T his makes th i s one 
of the few SMGs known at z > 4 (iCapak et al.l 1200 ' 
Schinnerer et al.l 120081: ICoppin et alj l2009t iDaddi et al 



2009allbf ) . In addition, a faint counterpart (r = 22.06 AB) 
is det ected in both the Sloan Digital Sky Survey (SDSS 
DR7, lYork et al.1 [2()00i and the United Kingdom In- 
frared Deep Sky Survey rUKIDSS. lLawrence et al.l[2007h 
that has a photom etric redshift of zi cns = 0.69 ± 0.13 
(|Smith et al.l l20l"l[) . Altogether, the evidence favors a 
scenario in which the background source is a SMG at 
high redshift that is being gravitationally lensed by an 
object at interme diate redshift, c onsistent w ith the lens- 
ing hypothesis of INegrello etatl ([20071 . [20101) . 

Recent observations by the Submillimeter Array 
(SMA) have shown an elongation along the southeast — 
northwest direction, possibly an indication of interest- 
ing morphological f eatures on scales smaller than « 
2" (|Cox et al.l 120111 ). In this paper, we present high- 
spatial resolution (0'.'6) SMA imaging at 880 /im, Keck 
Adaptive Optics (AO) Xs-band imaging, and Gem- 
ini Multi-Object Spectrograph-South (GMOS-S) optical 
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spectroscopy of this object and demonstrate their utility 
for constraining a detailed model of the lens-source sys- 
tem. We use the SMA, Keck, and Gemini data to probe 
the source size and magnification factor as well as the 
luminous plus dark matter mass of the lensing galaxies. 
The magnification factor is a critical parameter, since it 
is needed to understand intrinsic properties of the back- 
ground SMG such as its IR luminosity (Lir, integrated 
over 8-1000 fim) as well as molecular gas and dust masses 
(M gas and Mdust)- G15v2.779 is an example of a gravita- 
tionally lensed system discovered in H- ATLAS that per- 
mits the simultaneous study of obscured star-formation 
at high redshift as well as the nature of light and dark 
matter in galaxies at intermediate redshift. 

When the H-ATLAS catalogs are complete, we expect 
to have ~ 300 candidate lensed systems. In addition, 
systems identified from th e Herschel Multi-t iered Extra- 
galactic Survey fHerMES: lQliver et al.ll2010l ) could bring 
this total to 500 such objects discovered by Herschel. 
Recent efforts using near-IR imaging to push to fainter 
sub-mm flux densities and grow the list of candidates 
up to ~ 2000 appear promising (jGonzalez-Nuevo et al.l 
2012). Using similar selection techniques, the South 
Pole Telescope (SPT) has identified a sample of « 40 
candi date lensed SMGs within the initial 87 deg 2 sur- 
vey (jVieira et al.l 120101 ). The final survey area will 
cover 2000 deg 2 and is expected to provide a sample of 
~ 1000 lensed SMGs. Due to the selection at 1.4mm, 
the SPT sample of lensed SMGs will be complemen- 
tary to the Herschel sample in the sense that it will 
be biased towards higher redshift or cooler dust tem- 
peratures. Together, both the Herschel and SPT sam- 
ples offer the opportunity to expand upon the legacy 
of work on strongly lensed galaxies over a similar red- 
shift range undertaken as part of the Center for As- 
troph ysics Arizona Space Telescope Lens Survey (CAS- 
TLeS; lMuhoz et al.lll998h the Cosmic Lens All-Sky Sur- 
vey (CLASS: iMvers et al.lf200l iBrowne et al.ll2003f) . and 
the Jodrell Ba nk Very Large Array gravitational lens 
survey C7VAS: iKmg fc Browne) llOM) by increasing the 
sample size of such systems by 1-2 orders of magni- 
tude. In comparison to lensed systems selected via opti- 
cal spectroscopy (e.g., the Sloan Lens Advanced Camera 
for Surveys or SLACS and the Baryon Acoustic Oscil- 
lation Survey Emis s ion-Line Lensing Survey or BELLS; 
iBolton et all 120081: iBrownstein et al.l l2012h . the (sub- 
)mm selection is highly complementary in that it identi- 
fies (1) lensed galaxies that are both more luminous and 
at higher redshift; and (2) lensing galaxies that span a 
wider range in optical brightness (in particular, they do 
not need to be bright enough to be detected in SDSS-III 
spectroscopy). 

Throughout this paper we assume 
H =71 km s" 1 Mpc" 1 , Q m = 0.27, and tt x = 0.73. At 
z = 4.243, this results in a scaling of 6.9 kpc arcsec -1 . 

2. OBSERVATIONS 

2.1. SMA Data 

SMA imaging data of G15v2.779 was initially obtained 
as a short observation conducted in the compact array 
configuration on 2 010 June 16 (tj n t = 3 hrs on-source in- 
tegration time, see ICox et al.|[201lD . These data yielded 
a robust detection of the target and provided a total flux 



density measurement at 880 /jm of Fsso/im = (90±2) mJy. 
The data also pinpointed the location of the Herschel 
source: RA 14:24:13.98, Dec +02:23:03.45 (J2000.0) and 
with a beamsize of 2'.'4 x l'/3 hinted at an elongation in 
the nort hwest — southeast direction (—32 ±4 degrees east 
of north; ICox et al.ll20Tl . 

Subsequent data in the very extended array configura- 
tion were obtained on 2011 January 4 and 6 (tint = 3 hrs 
in total, max baseline length of 475 m). Extended array 
configuration data were obtained to improve uv coverage 
on 2011 January 28 (t lnt = 2 hrs, max baseline length of 
226 m). Atmospheric opacity was low (T225GHZ < 0.08) 
and phase stability was good (phase errors less than 30°). 

We optimized the SMA single-polarization 345 GHz 
receivers for continuum detection by tuning the pri- 
mary local oscillator to 339.58 GHz and an intermedi- 
ate frequency coverage of 4-8 GHz, providing a total of 
8 GHz bandwidth (considering both sidebands). The ob- 
servations did not cover the nearest CO emission line, 
CO( J = 16 — 15) at about 351 GHz since that frequency 
is near the high end of the range of the SMA 345 GHz 
receivers and would have compromised the sensitivity of 
our continuum observations. 

We used the Interactive Data Language (IDL) MIR 
package to calibrate the uv visibilities. The blazar 3C279 
was used as the primary bandpass calibrator while Titan 
was used as the absolute flux calibrator. The nearby 
quasars 1337—129 (F S8 o^m — 1-7 Jy, 19 degrees from 
target) and 1512—090 (F S8 o^ m = 1.5 Jy, 17 degrees from 
target) were used for phase gain calibration. The quasar 
1458+042 (^880^111 = 0-15 Jy ; 9 degrees from target) was 
observed to provide an independent check of the reliabil- 
ity of the calibration, particularly phase transfer. 

For imaging, we used the Multichannel Image Recon- 
struction, Ima ge Analysis, and Display (MIRIAD) soft- 
ware package (jSault et al.l 119951 ) with natural weighting 
for maximum sensitivity. This resulted in an elliptical 
gaussian beam with a full- width half-maximum (FWHM) 
of 0'.'69 x 0'.'50 and position angle of 62.4 degrees east of 
north. 

Figure [1] shows the SMA image of this source (com- 
bining compact, extended and very extended array 
configurations) in red contours. The sub-mm emis- 
sion is resolved into two bright emission regions to 
the southeast and northwest of the lensing galax- 
ies. These emission regions have peak intensities of 
Sgso^m = (21 ± 1-0) mJy beam -1 and S , 88 o M m 
: 1 1) ± 1.0) mJy beam -1 , respectively (quoted la un- 
certainties do not include the estimated absolute cal- 
ibration uncertainty of 10%). The peaks are located 
at positions of (RA=14:24:14.006, Dec=+02:23:02.81) 
and (RA=14:24:13.938, Dec=+02:23:04.40), respectively 
(the uncertainty in the relative position of these peaks is 
w 0'.'05). In addition, there is a background of complex 
substructure that, together with the two bright emis- 
sion regions, sums to a total flux density of Fg,sop, m = 
90 ± 1.0 mJy and is likely to be produced by gravita- 
tional lensing. 

2.2. Keck Adaptive Optics Data 

We obtained a 1920 sec -Ks-band image of G15v2.779 

on 2011 April 13 (UT) as part of program ID 

C213N2L (PI: H. Fu; e.g., iFu et al.l l20ll with the 
Keck II laser guide-star adaptive-optics system (LGSAO; 
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Fig. 1. — Multi-wavelength imaging of G15v2.779 centered on 
RA 14:24:13.975, Dec +02:23:03.60 (the 880Mm emission centroid) . 
Red contours highlight SMA 880 (im dust continuum emission from 
the lensed SMG at z = 4.243 (drawn at -2, 2, 4, 6, ... times the 
lrj rms level, where a = 1.0 mjy). The size of the SMA beam 
(FWHM 0'.'69 X 0'.'50) is shown by the black hatched ellipse at the 
bottom left corner of the plot. The greyscale background shows 
Keck AO Kg-band imaging which has resolved the lens into two 
early type galaxies. 



Fig. 2.— GALFIT modeling of G15v2.779. Left: Keck AO 
i\S-band image (as in Figure [TJ. Middle: Best-fit GALFIT model 
(assuming n s = 4 for both lensing galaxies). Bight: Residuals 
obtained after subtracting the best-fit model from the Keck data. 
Both morphologies are consistent with early- type galaxies. 

while the two foreground lensing galaxies are detected at 
high significance. The AO imaging indicates the sec- 
ondary lens galaxy is located 0'.'025 east and 0'.'327 south 
of the primary. A Galfit decomposition of the two 
sources seen in the Keck image into Sersic components 
indicates that a de Vaucouleurs profile is appropriate for 
both the northern and southern galaxies. In fact, the 
best-fit models have n s > 4, but this is probably due 
to faint, large-scale fluctuations in the background sky 
level - hereafter, we assume n s = 4 for both galaxies 
for simplicity. Figure [H shows the best-fit Galfit model 
and the residuals after subtracting the model from the 
data. Table [T] contains the best-fit parameters from the 
Galfit modeling and their la uncertainties. These are 
underestimates of the true errors as they do not account 
for degeneracies between the parameters. The two lens- 
ing galaxies are highly compact, with effective radii of 
r lonsl = 0.84 ±0.01 kpc and r lcns2 = 0.91 ±0.03 kpc. The 
northern galaxy has K§ ~ 18.22 ± 0.17 (Vega mag) and 
the southern galaxy has Kg = 19.01±0.17 (Vega mag). 

The background source is not detected inside a 1'.'5 ra- 
dius, implying a 5er limit of K$ > 20.14 (Vega mag). This 
corresponds to a flux density 5a limit of Fk s < 5.8 /iJy. 
This level of faintness in the n ear-IR is frequently an in- 
dication of z > 4 sy stems (e.g. iDannerbauer et all 120021 : 
IDaddi et al.ll2009al ). If differential magnification effects 
are insignificant (i.e., fiK s ~ A'sso/mi), then this implies 
a 500/im to Kg-band flux density ratio of w 35000. 
At z = 4.243, this corresponds roughly to rest-frame 
100/im to [/-band. An interesting comparison example 
is Mrk 231, a heavily obscured ULIRG in the local uni- 
verse that has a 100/im to [/-band flux density ratio of 
w 10000. This value is about a factor of 3-4 lower than 
G15v2.779, suggesting that the obscuration in G15v2.779 
is extreme. In section 14. 1[ we will compare the obscura- 
tion in G15v2.779 with other z ~ 4 SMGs. 

2.3. Gemini GMOS-S Optical Spectroscopy 

Long-slit spectroscopic observations of G15v2.779 were 
taken using the Gemini GMOS-South instrument on the 
night of 2011 March 6, under photometric conditions as 
part of programme GS-2011A-Q-57 (PI: D. J. B. Smith). 
Two observations of 1800 seconds each were made 
through a 2" slit, using the R400 grating and the OG515 
blocking filter, with dithering both in the wavelength di- 
rection and along the slit to minimize the effects of bad 
columns and gaps between the GMOS-S chips. The cen- 
tral wavelengths for the two observations were 630 and 

24 http://www2.keck.hawaii.edu/inst/nirc2/Manual/ObserversManual0-ggihm, and flat field observations were interspersed be- 

25 http://www2.keck.hawaii.cdu/inst/nirc2/forReDoc/ P ost.observing/ 1 d # gg fi p fe h e observations at each wavelength setting, as rec- 



iWizinowich et al.l 120061 ). An R = 17.9 magnitude star 
26" south of G15v2.779 served as the tip-tilt reference 
star. The expected Strehl ratio is about 0.17 at the 
source position. We used the second generation near- 
infrared camera (NIRCz£3) at 0'.'04 pixel -1 scale and ex- 
ecuted a nine-point dithering pattern with 3" dithering 
steps. Three 80 sec exposures were obtained at each 
dithering position. The natural seeing was about Of! 7. 

We used our own IDL program to reduce the images. 
After dark subtraction and flat-fielding, sky background 
and object masks were updated iteratively. For each 
frame, after subtracting a scaled median sky, the resid- 
ual background was removed with B-spline models. In 
the last iteration, we discarded the three frames of the 
poorest image quality and corrected the NIRC2 geomet- 
ric distortion using the solution of P. B. CamerorF^I be- 
fore combining the aligned frames. The resolution of the 
final image is O'.'l in FWHM, as measured from the stel- 
lar source 4'.' 5 northwest of G15v2.779. Astrometry was 
determined from SDSS photometric sources inside the 
40" field of view and carries la uncertainties in an ab- 
solute sense of « 0'.'4 (in fact, we show in section [3] that 
the tightest constraints on the astrometry are derived 
directly from the lens modeling). The flux scale in the 
image was normalized such that the total flux of the two 
lensing galaxies matches that seen in the UKIDSS K- 
band data, where the two galaxies have a total magni- 
tude of K = 17.89 ± 0.17. 

The grayscale of Figured] shows the Keck AO .Kg-band 
image of G15v2.779. The backgound SMG is undetected, 



A Lens Model for a SMG At z = 4.243 



5 



TABLE 1 
GALFIT Lens Modeling Results 3 - 







Lens 1 


Lens 2 


m Ks 


(Vega) 


18.22 ±0.01 


18.98 ±0.02 


^lcns 


b 


4 


4 


^lcns 


(kpc) 


0.84 ±0.01 


0.91 ± 0.03 


e lcns 




0.29 ±0.01 


0.06 ± 0.02 


</>lens 


(dcg) 


62 ±1 


-47 ±14 



a Uncertainties do not reflect degenera- 
cies between the parameters or absolute 
calibration uncertainty in m,K s ■ 
b Sersic indices are fixed to be ni ens = 4 
(see text for details). 



ommended by the Gemini observatory. CuAr arc lamp 
exposures were taken for the purposes of wavelength cal- 
ibration, using the same instrumental setup as for the 
science exposures, and the spectral resolution obtained 
was «6.0A. A position angle of 330° East of North was 
chosen to place the slit along the major axis of the ex- 
tension seen in the SMA and PdBI observations, and the 
CCD was binned by 4 pixels in both the spectral and 
spatial directions. The long-slit data were reduced using 
the IRAF Gemini GMOS reduction routines, and follow- 
ing the standard GMOS-S reduction steps in the exam- 
ple taken from the Gemini observatory webpages. Since 
the primary aim of these observations was to obtain a 
spectroscopic redshift and measure of the Balmer/4000 
A break for the lensing source (or sources), flux calibra- 
tion was not performed. 

Figure [3] shows the optical spectrum of G15v2.779 ob- 
tained with Gemini GMOS-S. Because the slit width is 
2", both lensing galaxies are included in this spectrum. 
A strong break at observed-frame 6350 A is obvious in 
the spectrum. Two narrow emission lines on either side 
of this break coincide with atmospheric features and are 
likely not produced by astronomical sources. The break 
in the spectrum likely corresponds to a Balmer/4000 A 
break and suggests that both lensing galaxies are at 
z = 0.595 ± 0.005 (the redshift and error were esti- 
mated using a by-eye comparison of the observed spec- 
trum with a synthesized simple stell ar population with 
solar met allicity and an age of 5 Gvn IBruzual fc Charlotl 
2003). We note that without spatially resolved spec- 
troscopy it is not possible to confirm that both lensing 
galaxies are located at the same redshift. The ratio of 
the flux density longward and shortward of the break 
is known as the D n (4000) value and is a measure of 
the average age of the stars within the gal axies. We use 
the definition from Kauffmann ct al. (2003) and measure 
D n (4000) = 2.0 ±0.2 (we adjust the wavelength window 
longward of the break to avoid the portion of the spec- 
trum affected by atmospheric emission lines). This value 
is typical of galaxies dominated by old s tellar populations 
(5-10 Gyr, depending on metallicity; Kauffman n et al.l 
120031) . 

3. A DETAILED LENS MODEL 

The combination of good sensitivity and high spa- 
tial resolution provided by the SMA and Keck AO data 
(probing the emission from the lensed and lensing galax- 
ies, respectively) permits a detailed study of the param- 
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Fig. 3.— Gemini GMOS-S optical spectrum of G15v2.779. 
Narrow cyan hatched lines indicate regions of the spectrum cor- 
rupted by either chip gaps in the CCD or atmospheric spectral 
features. Broad red-hatched regions indicate the portions of the 
spectrum used to compute D n (4000) and are selected to avoid the 
cyan regions. The background source is hidden by dust at these 
wavelengths, but both foreground lens galaxies lie within the 2" 
GMOS-S slit. The spectrum shows a strong, abrupt break at an 
observed-frame wavelength of A Q b s = 6350 A, likely correspond- 
ing to the Balmer/4000 A break and implying a lens redshift of 
2 lens = 0.595 ± 0.005. The strength of the break implies the lens- 
ing galaxies are dominated by old stars (5-10 Gyr, depending on 
metallicity; Kauffmann ct al. 20031). 

eters of the lensing model that describe this system. The 
lensed emission is detected only in the SMA data. Since 
the SMA is an interferometer, the surface brightness map 
of the lensed emission is obtained with incomplete uv 
coverage, implying that surface brightness is not neces- 
sarily conserved and that the pixel-to-pixel errors in the 
map are correlated. Furthermore, the lensed emission 
observed by the SMA comprises multiple, resolved com- 
ponents. For these reasons, it is important to compare 
model and data in the visibility plane rather than the 
image plane. 

We make use of t he publicly available Gravlens soft- 
ware (jKeetonl 12001) to map emission in the source plane 
to the image plane for a given lensing mass distribution. 
Using this software, a model of the lensed emission is 
constructed based on the density profile of the lens - as- 
sumed here to be two singular isothermal ellipsoids (one 
for each galaxy seen in the Keck AO imaging), the mor- 
phology of the source - taken here to be a Sersic profile, 
and the position of the source relative to the lens. Al- 
though the Sersic profile represents a crude simplification 
of the true background source morphology, we use it here 
because it permits a test of the variation in the lensing 
properties (e.g., magnification of the background source) 
as a function of variation in the nature of the source (e.g., 
half-light radius). Later in this section, we discuss how 
our best-fit parameters change when a second source is 
added to the source plane. 

The Gravlens lens model (with a single source in the 
source plane) contains 15 free parameters: the position of 
the source relative to the SMA 880/im emission centroid 
(Aa s and A<5 S ), the intrinsic flux of the source (F s ), the 
half-light radius of the source (r s ), the Sersic index of the 
source (n s ), the ellipticity (defined as srdi) anc j position 
angle of the source (e s and <fi s ), the position of the pri- 
mary lens relative to the SMA emission centroid (Aai ens i 
and A^ionsi), the mass of the primary lens (parameter- 
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ized in terms of the angular Einstein radius, #ei), the 
ellipticity and position angle of the primary lens (ei en si 
and Eiensi), and the Einstein radius, ellipticity, and posi- 
tion angle of the secondary lens (6*E2, £icns2, and </>i en s2)- 

The 1-tr absolute astrometric accuracy between the 
SMA and Keck images is 0'.'4, so in our modeling efforts 
we allow the position of the lens to vary by as much as 
0'.'8 in both RA and Dec (i.e., 2a in each direction). In 
fact, the constraints from the lens model are significantly 
tighter than 0'.'4, as we show below. 

We make use of the Keck AO .Kg -band data and fix 
the position of the secondary lens relative to the pri- 
mary at Aaiens2 = 0'.'025 and A£i ens2 = — 0'.'327 (rota- 
tional astrometric uncertainties between the SMA and 
Keck AO images are sub-dominant to translational un- 
certainties). Furthermore, we use our GALFIT results 
and constrain the ellipticity and position angle of the 
lenses to be within 3a of the best-fit Galfit values (as 
given in Table [TJ We tested models in which the mass 
of each foreground lensing galaxy was allowed to vary as 
well as models in which the mass of the secondary lens 
was fixed to be equal to one-half that of the primary lens. 
This latter choice assumes that both lens galaxies are lo- 
cated at the same redshift zi ons and is supported by the 
ground-based Gemini GMOS-S spectrum of G15v2.779 
(see section |2~B")) . High-spatial resolution optical or near- 
IR spectroscopy are needed to resolve the two lensing 
galaxies and provide a definitive test of the validity of 
our assumptions. 

For a given set of model parameters, Gravlens gen- 
erates a surface brightness map of the lensed emission. 
This surface brightness map can then be used as input to 
MIRIAD's uvmodel task, which computes the Fourier 
transform of the image and samples the resulting visi- 
bilities in a way that matches the sampling of the ac- 
tual observed SMA visibility dataset (Vsma) to produce 
a "simulated visibility" dataset (Vmodei)- The quality 
of fit for a given set of model parameters is determined 
from the chi-squared value (x 2 ) according to the follow- 
ing equation: 

2 _ (V S ma(u,v) - V mo dc\{u, v)) 2 



where a(u, v) is the la uncertainty level for each visibility 
and is determined from the system temperatures (this 
corresponds to a natural weighting scheme). Because 
the measured visibilities are complex, we compute both 
Xreai an d Xhnag an d measure the total x 2 as the sum of 
the real and imaginary components. 

We employ a Markov chain Monte Carlo (MCMC) 
technique to sample the posterior probability density 
function (PDF) of our model parameters. In particu lar, 
we use the emcee code (|Foreman-Mackev et al.ll2012h to 
implement the MCMC analysis. Th e algorithm adopted 
by em cee was originally presented in lGoodman fc Weard 
(2010) and uses an affine-invariant ensemble sampler to 
obtain significant performance advantages over standard 
MCMC sampling methods. 

Here, we summarize the behavior of the ensemble sam- 
pler technique. The available parameter space is searched 
using a set of iV wa ikers walkers. During iterations of the 
MCMC, each walker selects another walker from the en- 



semble and identifies a new position in parameter space 
based on the positions of both walkers (this is known as 
a "stretch move"). Once a new position has been found, 
the posterior PDF is computed and compared to that of 
the previous position. New positions with higher prob- 
ability (i.e., lower x 2 ) are always accepted; those with 
lower probability are sometimes accepted. After a suffi- 
cient number of iterations (A^ter), the ensemble of walk- 
ers samples the parameter space in a way that reflects 
the posterior PDF. The mean and variance of each pa- 
rameter can then be measured directly from the history 
of walker positions. We employ a "burn-in" phase with 
-^walkers = 250 and A^ter = 200 (i.e., 50,000 samplings of 
the posterior PDF) that is used to identify the best-fit 
model position. This position is then used to initialize 
the "final" phase with N wa ik eTS = 250 and Ai te r = 40 
(i.e., 10,000 samplings of the posterior PDF). To ensure 
that the posterior PDF was sampled with a sufficient 
number of iterations, we computed the autocorrelation 
time for each parameter in a given ensemble of walkers 
and found that it was of order unity for each parameter. 
This implies that we have 10,000 independent samplings 
of the posterior PDF, more than enough to obtain a ro- 
bust measurement of the mean and uncertainty on each 
parameter of the model. 

During each MCMC iteration, we also measure the 
magnification factor /i (we follow the nomenclature in the 
SMG literature here and use \i to refer to the total mag- 
nification obtained by summing over all individual lensed 
components) using the following method. First, we take 
the unlensed, intrinsic source model and measure the to- 
tal flux density (.Fin) within an elliptical aperture (^m) 
centered on the source with ellipticity and position angle 
equal to that of the source model and with a semi-major 
axis length of 5". Next, we take the lensed image of the 
best-fit model and measure the total flux density (F out ) 
within the aperture A ont , where A out is determined by 
using Gravlens to map A- m in the source plane to A out 
in the image plane (using the lens parameters which cor- 
respond to the best-fit model) . The magnification is then 
computed simply as /i = F out /Fi n . The best-fit value and 
la uncertainty are drawn from the posterior PDF as with 
the other parameters of the model. 

The best-fit model (assuming a fixed mass ratio be- 
tween the primary and secondary lens) produced by 
Gravlens is shown in Figure |4] and demonstrates that 
many of the features present in the SMA imaging can be 
reproduced in detail by Gravlens. The panel on the 
left shows the SMA imaging overlaid on the inverted, 
deconvolved map of the best-fit model visibilities. The 
panel on the right shows the residual image obtained by 
inverting and deconvolving the residual visibilities (i.e., 
the cleaned map of the difference between the model and 
data visibilities). The model fits the two brightest com- 
ponents in the SMA image (the peaks to the southeast 
and northwest), but fails to reproduce fully the peaks in 
the map to the northeast and southwest. 

We measure the following parameters of interest from 
the model: fi = 4.1±0.2, E i = 0'.'57±0'.'01, 6» E2 = 0'.'40± 
O'.'Ol, Aa lcns i = -0'.'27±0'.'03, A<5 lonsl = 0'.'63±0'.'03 (this 
position is within la of the position indicated from the 
Keck AO astrometry), Aa s = 0'.'03±0"02, A5 S = tf.'10± 
0'.'02, n s = 2.9 ± 0.3, r s = 4.4 ± 0.5 kpc, e s = 0.27 ± 0.09, 
and 4> s = 77±12 degrees east of north. The best-fit model 
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TABLE 2 

Gravitational Lens Model Results 





Single-source 
0El/0E2 free 


Single-source 
#El/0E2 fixed 


Two-source 
#Ei/0E2 fixed 




0.04 ± 0.02 


0.03 ± 0.03 


—0.05 ± 0.02 


AA 1 (") 


no _i_ n no 


n i n _|_ n no 


n no -|- n m 

VJ . WC< _1_ VJ . \J' J 


y i \\tY\r\ 
1 si VP/ 


3.2 ± 0.5 


4 4 -j- 0.5 


3.9 ± 0.6 




2.9 ± 0.3 


2.9 ± 0.3 


2.8 ± 0.5 




0.24 ± 0.09 


0.27 ± 0.09 


0.38 ±0.09 


4>sl (deg) 


84 ± 16 


77 ± 12 


80 ± 10 


Aa s2 (") 


— 


— 


-0.2 ±0.1 


AS s2 (") 


— 


— 


0.53 ±0.09 


r s2 (kpc) 






2.2 ±0.8 








0.5 a 








0.3 ±0.1 


0s2 (deg) 






90 ±30 


Aai onsl (") 


-0.28 ±0.03 


-0.27 ± 0.03 


-0.32 ±0.04 


A<5 lcnsl (") 


0.65 ± 0.04 


0.63 ± 0.03 


0.53 ± 0.04 


Obi (") 


0.44 ± 0.09 


0.57 ± 0.01 


0.57 ± 0.02 


e E2 (") 


0.53 ± 0.09 




e E1 /V2 




4.2 ±0.3 


4.1 ±0.2 


3.7 ±0.3 


x 2 


217796.9 


217799.1 


217769.2 


^DOF 


144190 


144191 


144185 



a Sersic index of secondary source is fixed at n s2 = 0.5. See 
sect ion [3] for details. 



has x 2 — 217799.1 and 144191 degrees of freedom. These 
constraints on model parameters are reported in Table O 
Also reported in Table [2] are the results obtained when 
the mass of the secondary is allowed to vary relative to 
the primary. The uncertainties on the mass of each indi- 
vidual lensing galaxy are significantly larger in this case 
because our lens model constrains the sum of the masses 
of the two lensing galaxies. The constraint on the sum 
of the masses is #ei + #E2 = 0'.'97 ± 0'.'02 (consistent with 
the results from models where the mass ratio has been 
fixed). The Gemini GMOS-S optical spectrum provides 
evidence that the two lensing galaxies are both z = 0.595. 
Under that assumption, the secondary is significantly less 
luminous and hence likely to be significantly less mas- 
sive. Therefore, for the remainder of the paper we have 
assumed a 2:1 mass ratio between the primary and sec- 
ondary. Spatially resolved spectroscopy of the lensing 
galaxies is needed to prove the validity of this assump- 
tion. 

The residual map in Figure |4] shows that our single- 
source model fails to account for the secondary peaks 
in the map to the northeast and southwest as well as a 
faint ring of emission partially lined up with the tangen- 
tial critical curve. We investigated whether adding an 
additional source in the source plane near the tangen- 
tial caustic would improve the ability of the model to 
match these secondary peaks. The new model has six 
new parameters: the flux of the second source {F S 2), the 
position of the second source relative to the SMA emis- 
sion centroid (Aa S 2 and A5 S 2, constrained to be within 
0'.'2 of the tangential caustic), half-light radius of the sec- 
ond source (r S 2), the ellipticity of the second source (e S 2), 
and the position angle of the second source (4> S 2)- Initial 
tests of the two-source model showed that the fitting rou- 
tine struggled to identify the best-fit model consistently 
unless the Sersic index of the second source was fixed at 
n S 2 = 0.5 (i.e., a Gaussian profile), so that is what we 
have adopted here. 



Figure[5]shows the results obtained with the two-source 
model and Table [2] contains the best-fit model parame- 
ters and their la uncertainties. The best-fit model has 
X 2 = 217769.2 and 144185 degrees of freedom. The mag- 
nification factor is fi = 3.7±0.3, and the Einstein radii of 
the two lenses are E i = 0.57±0.01 and 9 E2 = 0.40±0.01, 
values that are very similar to those obtained with the 
single-source model. These are the parameters of great- 
est interest, so it is reassuring that they are relatively 
insensitive to the exact morphology of the background 
source. 

The position of the primary background source is also 
relatively robust between the two models we have tested, 
with (Aq s = 0'.'05 ± 0'.'02, A5 S = 0'.'09 ± 0'.'03) relative to 
the SMA emission centroid. The secondary source is less 
well-constrained, having a position of Aa S 2 — 0'.'2 ± ff.'l, 
A5 S 2 = 0'.'53 ± 0'.'09. This position places the secondary 
source near a caustic, implying that it has experienced 
a high degree of magnification and has a much lower in- 
trinsic luminosity than the primary source. The primary 
source maintains a broad profile intermediate between a 
disk and an elliptical, while the second source is compact 
(r S 2 = 2.2 ± 0.8 kpc). In both the single-source and two- 
source models, the best-fit residual maps show emission 
to the southwest in the SMA data that the model fails to 
reproduce. This could be an indication of the presence 
of a third source in the source plane, but an examination 
at that level is beyond the scope of this paper. Overall, 
the agreement between the single-source and two-source 
model results is encouraging. For the remainder of the 
paper, we use the results from our single-source model, 
since it provides nearly as good a fit to the data as the 
two-source model (A% 2 = 29.9) but requires six fewer 
free parameters. Our major conclusions are insensitive 
to whether the single-source or two-source model is used. 

It is interesting that the magnification factor we have 
measured is somewhat lower than might be expected 
based on the velocity dispersion and observed lumi- 
nosity of CO emission lines in G15v2.779. Studies of 
CO (J = 1 — 0) emission lines in both lensed and un- 
lensed SMGs have found a correlation between the in- 
trinsic (i.e., unlensed) CO(J =1 — 0) line luminosity and 
the FWHM of the emission line (Harris et al., submit- 
ted; Bothwell et al., in prep.). While CO(J =1-0) 
data are not yet available for G15v2.779, higher J lines 
have indicated FWHM= 700 km s" 1 and ^co(J=5-4) = 
3xlO n A'kms -1 pc 2 . According to the correlation iden- 
tified in the previous studies, this FWHM value would 
imply an (unlensed) line luminosity of ^co(j=i— o) = 
7 x 10 10 K kms -1 pc 2 . However, assuming a typical 
sub-thermal CO(J = 5 - 4) to CO(J =1-0) con- 
version factor (e.g., Harris et al., in submitted.), our 
measurement of /i would indicate ^co(J=i— o) = ^ x 
10 11 K kms -1 pc 2 . Possible explanations for this in- 
clude either an unusual CO (J = 5 - 4)/CO(J =1-0) 
ratio, a large intrinsic scatter in the luminosity line-width 
relations in these systems, or a difference in the size-scale 
of the CO (J =1 — 0) and the far-IR emitting regions, 
leading to /ico not being equal to /if ar -iR- 

4. THE NATURE OF G15V2.779 

In this section, we discuss the implications of our re- 
sults for the nature of G15v2.779. We begin with a focus 
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Fig. 4. — Comparison of best-fit Gravlens model with SMA data. Left: SMA imaging (red contours) overlaid on the inverted, 
deconvolved map of the best-fit model visibilities. Right: Residual image obtained by inverting and deconvolving the residual visibilities 
(i.e., the cleaned map of the difference between the model and data visibilities). Critical curves and caustics are traced by orange and cyan 
lines, respectively. Lens positions as derived from the Keck AO i-fg-band imaging are marked by plus signs. The peak flux position of the 
source is shown with a filled blue circle. Contours indicate -2, 2, 4, 6, etc. times the la rms noise level. This model has a magnification 
factor of n = 4.1 ± 0.2 and Einstein radii of E1 = 0'.'57 ± O'/Ol and E1 = 0'/40 ± O'.'Ol. 
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Fig. 5. — Same as FigureU but now with a two-source model (filled blue circles indicate the peak flux positions of the two galaxies in 
the source plane). This model has a magnification factor of fi = 3.7 ± 0.3 and Einstein radii of = 0'/57 ± 0'.'02 and 6^2 = 0"40 ± O'.'Ol. 
This suggests that the measurement of these important parameters is relatively insensitive to whether a single-source or two-source model 
is assumed. 



on the background source and end with the foreground 
lens properties. 

4.1. The Background Source: An optically-obscured 
SMG at z = 4.243 

The magnification factor of the background source is a 
parameter that is critical to the derivation of any intrinsic 
property of the source. The properties that we consider 
here are the total IR luminosity (£tr), the SFR, the pro- 



jected IR luminosity surface density (Sir), the dust mass 
(A^dust), and the molecular hydrogen gas mass (Mh 2 ). 
We comput e £ir as the integ ral under the full SED as 
reported in ICox et all ([201 ll ). divided by \i. Differen- 
tial magnification effects are an important consideration 
when computing Xir . Initial studies indicate that when 
the magnification factor at any given wavelength is mod- 
est (/j, < 5), then differential magnification effects are 
expected to have a minimal influence on the inferred IR 
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luminosity (Serjeant 2012|). The contribution from the 
lensing galaxies to Lir is likely to be minimal, given that 
they are both early- type galaxies at z = 0.59. We find 
an intrinsic IR luminosity of Lir = (2.1 ±0.2) x 10 13 L Q , 
a value that likely makes G15v2.779 one of the intrinsi- 
cally brightest known SMGs (with an intrinsic, un-lensed 
sub- mm flux density of i^so/nm (intrinsic) = 21 ± 2 mJy). 

Using the sta ndard conversion from Lir to SFR 
(|Kennicuttiri998t ). this corresponds to a SFR of 3000 ± 
300 Mq yr" 1 . This value assumes that all of the IR lumi- 
nosity originates from star-formation rather than AGN. 
Hence this value should be regarded as an upper limit on 
the true SFR, although it should be noted that measure- 
ments of the radio luminosity suggest a far-IR to radio 
flux rati o that is consist ent with starburst dominated 
galaxies (|Cox et al.ll2011[ ). 

The Sersic index of the background source (n s = 2.9 ± 
0.3) is intermediate between an exponential disk profile 
(n s = 1) and a De Vaucouleurs (n s — 4) profile. The 
source is significantly extended, having a half-light radius 
of 4.4±0.5 kpc. This implies a de-projected IR luminosity 
surface density (computed as Sir ~ 0.5LiR/Ah a if , where 
A hal{ = 7rr s 2 ) of Sir = (3.4 ± 0.9) x 10 11 Lg ) kpc" 2 . This 
range of values for Sir places G15v2.779 within the "nor- 
mal star-forming" class of galaxies at z > 2 and is signif- 
icantly below that of local ULIRGs (which have Sir val- 
uesjn the range 10 12 — 10 14 L Q kpc~ 2 ; Ruiopakar n et al.1 
2011). This may be a clue that unlike local ULIRGs, the 
intense, dust-enshrouded burst of star-formation that is 
occurring in G15v2.779 may not be driven by a major 
merger. 

To compute the apparent ( i.e., u n-corrected for lens- 
ing) dust mass, ICox et all ()2011f ) perform a single- 
temperature, optically thin, modified black-body fit to 
the far-IR and su b-mm SED of G15v2.779 (following 
iBeelen et al.l [2006) . These authors find a dust mass 
of M dust = 8.9 x 10V -1 M Q , using the best-fit dust 
temperature of Ldust = 38 K and a mass absorption 
coefficient of k = 0.4 cm 2 g _1 . If instead an opti- 
cally thick modified black-body is used to fit the data, 
the inferred dust mass decreases by nearly a factor of 2 
(jWeifi et al.l 120071 ). Our measurement of fi implies dust 
m asses M^, K t ~ 2 x 10 9 Mq (using the values reported 
in ICox et all 1201 lh . The uncertainty in the dust mass 
is dominated by systematic uncertainties related to the 
unknown optical depth at far-IR wavelengths and the 
unconstrained mass opacity coefficient and is at least 
a factor of a few. The measured dust mass is similar 
(though at the top end) to th ose estimated for other high- 
z SMGs (jSantini et al.ll2010h and a lso for the most m as- 
sive z = 0.5 galaxies in H- ATLAS (|Dunne et al.ll2011l ). 

To quanti fy the int e rstella r medium properties in 
this object, ICox et all (|2011l) use a spherical, single- 
component, larg e velocity gradient (LVG) model 
(jWeifi et all 120071 ) to fit simultaneously the CO( J = 
7 - 6), CO(J = 5-4), and CO(J = 4-3) emis- 
sion lines. Assuming a conversion factor of aco = 
0.8 Mq (K km s -1 pc 2 ) -1 to go from Lco(i-o) to 
M H2 , ICox et al.l (l20ll find gas masses of Mh 2 = 3.5 x 
10 11 /i -1 Mq. Using our measurement of /i, we find 
M Ha = (8 ± 4) x 10 10 Mq. Note that there is some ev- 
idence that aco may be higher for less dense sy stems 
at high-z relative to local ULIRGS (jFraver et al.l 120111: 



llvison et al.l 1201 lT) . If this is indeed the case, then the 
gas mass would be even larger. Given the large dust 
mass (even with the optically thick fit), this is not un- 
reasonable. 

This gas mass is a factor of « 2.5 lar ger than the typica l 
gas mass found in un- lensed SMGs (IGreve et al.l 120051) . 
2.5-9 times larger t han two lensed sou rces discovered in 
the H- ATLAS SDP graver et al.ll20rl . and a factor of^ 
2.5 la rger than a lensed source from HerMES (| Scott et al.l 
l20TlT ). G15v2.779 appears to be a very massive, highly 
star-forming galaxy at z = 4.243. 

Table [3] summarizes the properties of G15v2.779 and 
compares them with unlensed SMGs at z ~ 4. G15v2.779 
bears a close resemblence to GN20, having similar Lir, 
Mg &s , (Tgas, r a , Mdyn, and Sir values. However, there 
is an important difference: GN20 is much more lumi- 
nous in the rest-frame optical than G15v2.779. In terms 
of visual extinction (Ay), G15v2.779 is most similar to 
GN10, which has A v ~ 5 - 7.5. However, GN10 is less 
luminous and has only an upper limit on its source size. 
The fact that G15v2.779 is clearly extended on scales 
> 2 kpc and yet maintains a covering fraction near unity 
indicates an impressive quantity of dust, consistent with 
the dust mass measurements described above. 

One feature all of these z ~ 4 SMGs share in com- 
mon is Sir values that are 1-2 orders of magnitude lower 
than those of ULIRGs in the local Universe. This may 
be an indication that the physical mechanisms driving 
the prodigious luminosities in these systems may be dif- 
ferent from what occurs at z ~ 0. Theories attempting 
to explain the behavior of these systems must also ac- 
count for the short gas depletion timescales (possibly via 
strong inflow of ga s from the inter-galactic medium; e.g. 
iSchaye eTa l. 2010), as the estimated SFRs will consume 
all of the available gas within ~ 10 — 30 Myr in all of 
these galaxies. 

4.2. The Foreground Lenses: A Dry Merger at 
z = 0.595 

We use the standard equations from iSchneider et al.l 
(1992) to compute the mass of the lens galaxies Mi ens , 
finding Mi ens i = (7.4 ± 0.5) x 10 10 Mq (recall that we 
have assumed Mi cns i = 2 x Mi cnS 2). Independent mass 
estimates can be obtained using the corre lation between 
U-ban d luminosity and Mi cns given in Nc grello et al.l 
(2010). At zi cns = 0.59, the observed z-band corre- 
sponds almost exactly to the rest-frame U-band. The 
ground-based z-band mag nitude of the t wo lenses to- 
gether is z = 20.35 ± 0.40 (|Cox et al.ll2011l ). correspond- 
ing to a rest- frame U-band luminosity of vL u (V) = 
3.3 jz^'^ x 10 10 Lq. Assuming the rest-frame U-band lu- 
minosity ratio between the two lens galaxies is the same 
as that in the observed i^s-band, the Ly — Mi fins correla- 
tion observed from th e SLACS lenses (jBolton et al.ll2008l: 
iNegrello et al.l I2010D implies lens masses of Mi ens i = 
(4 ±2) x 10 10 Mq and M lcns2 = (2±1) x 10 10 M , consis- 
tent with our lens model estimates of Mi cns i and AL. onS 2. 
An alternative estimate of the lens masses is to use the 
U-band luminosities along with a mass-to-light ratio de- 
rived from synthesized stellar populations. The optical 
spectrum of G15v2.779 shows a strong Balmer/4000 A 
break (D„ (4000) = 2.0 ± 0.1), typical of galaxies domi- 
nated by old stellar populations (5-10 Gyr, depending on 
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TABLE 3 
Properties of z > 4 SMGs 





^IR 
(10 13 Lq) 


SFR 
(M Q yr" 1 ) 


/agn 


M 

11 1 gas 
(10 11 Mq) 


A^stars 

(10 11 Mq) 


A v 


& gas 

(km s- 1 ) 


(kpc) 


M dyn sin 2 i 
(10" Mq) 


Sir 

(10 12 Lq kpc" 2 ) 


G15v2.779 


2.1 ± 0.2 


3000 ± 300 


low fg 


0.8 ±0.4 




> 4 


800 ± 100 


4.4 ±0.5 


3±1 


0.34 ± 0.09 


GN10 a 


1.2 ±0.6 


1700 ± 800 


low s 


0.27 ±0.05 


1.0 ± 0.5 


5 - 7.5 


770 ± 200 


< 4 


< 2.5 


> 0.1 


GN20 b 


2.9 ± 1.6 


4300 ± 2000 


low g 


0.50 ±0.25 


2.3 ± 1.5 


< 2 


710 ± 120 


2± 1 


2.3 ± 1.5 


1.1 ± 0.5 


GN20.2a b 


1.6 ±1.0 


2300 ± 1100 


low s 


0.30 ±0.15 


0.5 ±0.3 


< 2 


1100 ±400 


< 4 


< 5 


> 0.1 


J1000+0234 c 


1.2 ±0.7 


1700 ± 900 


low g 


0.26 ±0.13 


0.30 ± 0.15 


1.4 ±0.5 


600 ± 200 


3.5 ±2.0 


1.3 ±0.7 


0.2 ±0.1 


J033229.4 d 


0.6 ±0.3 


900 ± 450 


low h 


0.16 ±0.03 


< 0.5 


1.5 ±0.5 


160 ± 65 


2± 1 


0.12 ±0.06 


0.2 ± 0.1 


AzTEC-3 e 


1.7 ±0.8 


1800 ± 900 


low g 


0.53 ±0.25 


0.10 ± 0.05 


< 2 


487 ± 58 


< 5 


> 1.4 


> 0.1 



a IDaddi et~aH d2009al> 

b Daddi ct al. (2009b); Carilli ct al. (2011) 
c Capak ct al. (2008); Schinncror ct al. (2008) 
d ICoppin etall f200g| [207ol) 
e IRiechers et al.l (I2010l 'l 
f 880/^111 size > 2 kpc 

g Radio to IR. luminosity ratio consistent with star-formation dominated galaxies 
b Weak NV, 24/im, and near-UV emission inconsistent with AGN 



metallicitv; Kauffmann et al.ll2003|) . A synthesized stel- 
lar population with an age of 5 Gyr, a Chabrier IMF, 
and no dust extinction has a mass-to-light ratio of 
Mst^/vLy (V) k, 2. Thus, the inferred stellar masses of 
the two lensing galaxies are M sta r,i ~ (4 ± 2) x 10 10 M Q 
and Mstar,i ~ (2±1) x 10 10 M Q . There is evidence at the 
1.5a level that the lens and mass estimates are greater 
than the stellar mass estimates. This may be reason- 
able considering the small stellar sizes (half-light radii of 
« 0.9 kpc) relative to the Einstein radius of each system 
(~ 3 kpc). 

Figure [5] shows that G15v2.779 has two of the faintest 
(in i-band) and lowest mass (at this redshift) lensing 
galaxies found in current surveys of gravitationally lensed 
systems (e.g., CASTLeS, CLASS, SLACS, and BELLS; 
Mufioz et al1ll998l: IMvers et al.|[200l IBolton et aT1[200l 
Brownstein et al.l |2012). This is an indication that wide- 
field surveys with Herschel will be useful for identify- 
ing lensing systems where the lensing galaxy is faint in 
the optical, either due to being low-mass or being very 
distant. This is generally true for source-selected lens- 
ing surveys (e.g., CLASS, H-ATLAS, and a portion of 
CASTLeS), whereas lens-selected lensing surveys (e.g., 
SLACS and BELLS) tend to be biased towards brighter 
and more massive foreground galaxies. 

Both G15v2.779 lens galaxies have early-type mor- 
phologies with small half-light radii (« 0.9 kpc). Such 
systems appear to be commonplace at z ~ 2, but become 
incre asingly rare at lower redshifts (e.g.. lDamjanov et al.1 
2009). However, given the small separation between the 
two lensing galaxies (« 2 kpc), it is likely that they are 
about to merge together (dissipationless mergers often 
do not have obvious signs of interac tion even at these 
separations; e.g., see lBell et al.l 120061 ). Simple virial ar- 
guments suggest that this process could lead to a dou- 
bling of the radius while the mass only increases by 50% 
(|Naab et al.l 120 09) . Such a result would make the size 
of the merged system c onsistent with si milarly massive 
galaxies at z « 0.5 (e.g. JOser et al.ll2012D as well as more 
typical of lensing galaxies found in CASTLES, CLASS, 
SLACS, and BELLS. Overall, these observations are con- 
sistent with the dissipationless ("dry") merging stage 



of the two - phase galaxy evolution scenario outlined in 
lOser et all (|2010( ). In this picture, the progenitors of 
massive galaxies undergo intense in-situ star-formation 
from z ~ 6 to z ~ 2 that leads to compact, elliptical 
galaxies with little or no reservoirs of gas for future star- 
formation. From z ~ 2 to the present-day, massive galax- 
ies undergo dry merging and increase their sizes such that 
they evolve onto the mass-size relation observed in local 
early-type galaxies. 

5. CONCLUSIONS 

We use high-spatial resolution imaging obtained with 
the SMA at 880/mi and Keck AO at if s -band to per- 
form a detailed gravitational lens modeling of G15v2.779, 
an SMG at z = 4.243 identified by Herschel in the H- 
ATLAS survey. We present a Gemini GMOS-S optical 
spectrum of G15v2.779 that suggests that the two fore- 
ground galaxies are at zi cns = 0.595 ± 0.005. This anal- 
ysis provides important measurements of the nature of 
both the background SMG and the foreground lenses. 
We summarize our findings below. 

We employ a visibility-plane lens modeling analysis 
and find a magnification factor of fi = 4.1 ± 0.2 for the 
background source. This measurement is significantly 
lower than what had been previously assumed for this 
source and indicates that not all of the brightest lens 
candidates identified by Herschel have high magnifica- 
tion factors. This value of /i implies an intrinsic infrared 
luminosity of Lir = 2.1 ± 0.2 x 10 13 L Q . 

The best-fit model for the background source favors 
radial profiles that are intermediate between exponential 
disks and De Vaucouleurs. The half-light radius of the 
background source is r s = 4.4 ± 0.5 kpc. This size mea- 
surement implies a de-projected IR luminosity surface 
density of E IR = (3.4 ±0.9) x 10 n L Q kpc" 2 . This num- 
ber is typical of z > 2 ULIRGs and HyLIRGs but 10-100 
times lower than ULIRGs in the local Universe. This 
may be an indication that the formation mechanism for 
this source could be different from z ~ ULIRGs, which 
are thought to arise from major mergers of gas-rich disk 
galaxies. Higher-spatial resolution data with improved 
sensitivity are needed to favor one of these models over 
the other. 
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Fig. 6. — Comparison of lens properties from BELLS 
(gray square s; Brownstein ct al. 2012), CASTLeS (red tri angles; 
IMurloz et all 119981 ). CLASS (green stars: Kers et al.ll2003f l, and 
G15v2.779 (black filled circles represent individual lens galaxies; 
open circle represents the sum of both lens galaxies). Top: Lens 
i-band magnitudes as a function of redshift. Note that G15v2.779 
is significantly fainter in i-band than any of the galaxies in BELLS. 
Bottom: Mass within the Einstein radius as a function of redshift. 
G15v2.779 is less massive than all BELLS galaxies at 2i cns > 0.5. 
This demonstrates that as a source-selected lens survey (similar 
to CLASS and part of CASTLeS), H- ATLAS will be less biased 
towards massive, bright lensing galaxies than lens-selected surveys 
like SLACS or BELLS. H- ATLAS and other Herschel wide-field 
surveys will help identify lens systems at intermediate redshifts 
where the galaxy in the foreground is particularly faint, either due 
to having a low mass or lieing at high redshift (or both, as appears 
to be the case for G15v2.779). 

Our measurement of /i, in conjunction with previous 
observations of CO emission lines and the far-IR SED, 
indicates a gas mass of Mh 2 ~ (8 ± 4) x 10 10 M© and 
a dust mass of Md us t ~ 2 x 10 9 M©. These values are 
factors of 2.5-9 times larger than other lensed galaxies 
studied to date but are comparable to other z ~ 4 SMGs. 
They indicate that G15v2.779 hosts a massive reservoir 
of molecular gas that is fueling a prodigious, but likely 
short-lived (~ 10 — 30 Myr) period of star-formation. 

The foreground lenses have Einstein radii of 6ei = 
0"57 ± O'.'Ol and 9 E2 = 0'.'40 ± O'.'Ol. These imply 
lens masses of A/i cns i = (7.4 ± 0.5) x 10 10 M and 
Micns2 = (3.7±0.3) x 10 10 M©. The lensing galaxies have 
sizes of ~ 0.9 kpc and lie at a redshift of z — 0.595±0.005. 
The gravitational potential from both galaxies may in- 
clude a significant contribution from dark matter. They 



are separated by 2 kpc and will likely merge into a sin- 
gle early-type galaxy with a larger size that will make 
the resultant system consistent with other early-types at 
z~0.6. 

Together, the SMA, Keck, and Gemini data have es- 
tablished that G15v2.779 is a SMG at z = 4.243 mod- 
estly lensed by a pair of early- type galaxies at z — 0.595. 
Our results highlight the bounty of information that can 
be obtained via a multi-wavelength approach to study- 
ing strongly lensed SMGs at high redshift. More sensi- 
tive and higher-spatial resolution imaging of the lensed 
emission is needed to improve the constraints on the pa- 
rameters of the gravitational lensing model and test com- 
peting models for the powering mechanism in this source 
(e.g., major merger vs. secular processes). This will be- 
come feasible in the near future when baseline lengths of 
w 1 km become available with ALMA. HST spectroscopy 
is needed to confirm that the lensing galaxies both lie at 
z = 0.59, to measure their velocity dispersions, and to 
improve their stellar mass estimates. 

The results described in this paper are based on ob- 
servations obtained with Herschel, an ESA space obser- 
vatory with science instruments provided by European- 
led Principal Investigator consortia and with important 
participation from NASA. The Herschel- ATLAS is a 
project with Herschel, which is an ESA space obser- 
vatory with science instruments provided by European- 
led Principal Investigator consortia and with impor- 
tant participation from NASA. The H- ATLAS website is 
http://www.h-atlas.org/. US participants in H- ATLAS 
acknowledge support from NASA through a contract 
from JPL. RSB acknowledges support from the SMA 
Fellowship program. HF, AC, JLW and SK acknowledge 
support from NSF CAREER AST-0645427. We thank 
the referee for a thorough review of the manuscript which 
resulted in a stronger paper overall. 

The ground-based follow-up observations were ob- 
tained at the SMA, at the W. M. Keck Observatory, and 
at the Gemini South Observatory. The SMA is a joint 
project between the Smithsonian Astrophysical Obser- 
vatory and the Academia Sinica Institute of Astronomy 
and Astrophysics and is funded by the Smithsonian In- 
stitution and the Academia Sinica. The authors wish to 
recognize and acknowledge the very significant cultural 
role and reverence that the summit of Mauna Kea has 
always had within the indigenous Hawaiian community. 
We are most fortunate to have the opportunity to con- 
duct observations from this mountain. 

Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Coun- 
cil (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia, Tecnologia e 
Inovacao (Brazil) and Ministerio de Ciencia, Tecnologia 
e Innovation Productiva (Argentina). Facilities: SMA, 
Keck, Gemini-South. 
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